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ABSTRACT Conditions sufficient for a quantum dynamical semigroup (QDS) to be unital are proved for a class 
of problems in quantum optics with Hamiltonians which are self-adjoint polynomials of any finite order in creation 
and annihilation operators. The order of the Hamiltonian may be higher than the order of completely positive part 
of the formal generator of a QDS. 

The unital property of a minimal quantum dynamical semigroup implies the uniqueness of the solution of the 
corresponding Markov master equation in the class of quantum dynamical semigroups and, in the corresponding 
representation, it ensures preservation of the trace or unit operator. We recall that only in the unital case the 
formal generator of MME determines uniquely the corresponding QDS. 
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Numerical experiments with the remarkable quantum trajectories algorithm [Q] (Schack & Brun '97), which solves 
the Markov master equation (MME) by the Monte-Carlo method, show that the numerical code forces a solution to 
be unital even in cases, when the exact minimal solution does not preserves the trace of initial state or unit initial 
operator, and the corresponding Poisson process explodes on finite time intervals. Unfortunately, the mathematical 
theory of MME still does not work for many basic equations in quantum physics. In this paper we extend the 
approach developed in ^ to MMEs in quantum optics. 

The formal generator £,{■) of a quantum dynamical semigroup on von Neumann algebra B{Ti.) of all bounded 
operators in a separable Hilbert space Tl is called regular if it defines the semigroup in an unambiguous way. 
Following 1^, we assume that the coefficients of the generator £{■) are a densely defined symmetric Hamiltonian 
operator H and a completely positive map "!'(•). The structure of the generator £,{■) is similar to the structure 
of the classical Kolmogorov-Feller equation, and similarly to the classical case, under rather general assumptions, 
MME has the minimal solution called a minimal quantum dynamical semigroup (QDS). Moreover, if the minimal 
QDS is unital, it is the unique solution of the corresponding MME; this implies regularity of the formal generator 
£(•). 

In the present paper we suggest a new test for regularity of MME with the formal generator £,{■). The idea consists 
in a suitable choice of some A-pair consisting of a "reference" operator A and "interaction" part Hint = A + Hs.a., 
which generates the interaction representation for the Hamiltonian H — Hs + Hs.a. = {Hs — A) + (A + Hs.a.), 
where Hs and Hs.a. are symmetric and self-adjoint components of H. 

We call the pair (A, Hs.a.) a A-pair for the generator £{■), if the positive self-adjoint A is a reference operator for 
the reduced generator Co{-), i.e. 

jCoiX) = <1>(X) - G;X - XGo, Go = 1^1) + i{Hs - A), 



£o(A) < cA, I < $(/) < A, 



and Hs.a. and Hs are such that 



-f^Hint < Hs.a < fiHint, < Hint < f^, 

< ciA, ^{Htnt) < C2A, 



HHs 

for some c, fj,, u > and e G (0, 1). Under these algebraic assumptions, together with some domain and continuity 
conditions which will be specified in Sections 1 and 2, we prove that the minimal dynamical semigroup with the 
formal generator £(•) is unital. 

The paper is organized as follows. In Section 1, we discuss the construction of a minimal solution for the Markov 
master equation with time-dependent coefficients and prove new rather general conditions sufficient for the con- 
servativity of a minimal solution. In Section 2, we discuss the properties of the interaction representation of MME 
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and introduce definitions of a reference family and A-pair. In Section 3, we describe the properties of generators 
of MME in quantum optics, and present basic examples of violation of the unital property. In Section 4 several 
physical examples of nontrivial MME are considered which in the most cases are regular. 

The main conclusion of the paper is rather simple: for the class of CP-maps $(■) in quantum optics, the minimal 
QDS is unital if the Hamiltonian H is a, s.a. operator. 



1 Conditions sufficient for conservativity 

Consider the Markov master equation with time-dependent generator 

d 



^^P.AB)=Ct{PrAB)), PrAB)\r=t^B, (1.1) 

Ct{B) = ^t(B) - GIX - XGt, Gt = i^tW + iHt. 
and assume that there exists strongly continuous contractive evolution system Ws,t £ S(7i) {t > s) such that 

Ws,rWr,t = Ws,t, i-Ws,t = -Ws,tGt, ^W,,t = GsW,,t. 

at as 

For simplicity, we assume that all the generators Gt have a joint core Dn C dom G't for some A*' > 1 such that 
the following preliminary domain assumptions are fulfilled 

V(/s,t : ©AT ^ ©AT C dom $.,((/)[■], G:+G. on Oat, (1.2) 

and that there exists some positive self-adjoint operator A : A > $«(/) > / such that the preliminary continuity 
conditions are fulfilled: 

the family of CP-maps A~^^^$t{-)^~^''^ is bounded, normal and ultraweakly continuous in t, and on the other 
hand, for any ip G T>m, the family of vectors ^{s,t) = X''^'^Ws,tip belongs to LfflM.+ jT-L) in variable s and is 
norm- continuous in variable t. 

Under these domain and continuity conditions, for any bounded strongly continuous family of operators Xa, the 
family of quadratic forms 

^s{Xs)[Ws,t^p] VV^eDiv (1.3) 

belongs to Li°'^(R+) in variable s and is continuous in t. 
Hence the sequence of CP-maps 

P^'hB) = W:,tBWr,t VrAB), 

Pi"t^'\B) VrAB) + J' dsVs,t^sVr,sPthB) 

is well-defined as a sequence of bounded operators corresponding to the sequence of densely defined and uniformly 
bounded quadratic forms. Indeed, K-,i(/) < I, and the identity 

VrAI)=I- dsVs,t^s{I) ^ I - J' dsVr,s$s{I) (1.4) 

readily shows that the sequence Pi"\B) is uniformly bounded: 

ll^'i"'(B)||<||i3||. 

Moreover, it increases monotonically if B G B+(Tt), and defines the least upper bound: 

-Pm"(-B) = l.u.b. Pi"' (B). 

This construction is analogous to the construction of the minimal solution for the Markov master equation with 
constant coefficients [^-[^]; it was discussed in details in 
The identity (1.4) implies that 

Pi]}{I) = VrAI) + ds Vs,t^sVrAI) 

= VrAI)+ I dsK,t$s(/)- / dSlKi,t$si / ' dS2Vs„s,^sM) 



/t rsi 
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Similarly, by using sequentially the identity (1.4), we obtain 

A'"'(r,t)=^^ dSlVs^t^s,--- J^"'' dSnVs^,s^_,^s„{I). (1.5) 

Hence, P^"j^{I) — > I strongly as n ^ oo, if and only if A„(s,t) — > weakly. 

To prove a condition sufficient for the minimal solution of the Markov master equation (1.1) to be unital, let us 
consider an estimate for the integral of the operator P^"'<1>t(/): 

ds Pi"> $,(/)= ^ dsVs^t^s{I) + ... 
+ I ds I dSl Ki,t$sr ■ ■ / dSnVs„,s„_i^s„Vs,s„^s{I). 



In the last multiple integral, the variables Sk take greater values then s, i.e. t < s < s„ < ■ ■ ■ < si. Hence by 
changing the order of integration, we have 

J' ds Pi"'$.(^) = J' ds Vs,t^s{I) + ■■■ 

+ [ dSlKi,t$si--- / " ' rfSn V;„,s„_i<&s„ / " dsVs,s„^s{I) (1.6) 



The last integral in (1.6) can be rewritten as the last integral in (1.5) in notation s — > s„+i. By comparing Eqs. (1.5) 
and (1.6) and passing to the least upper bound in n, we obtain the following important equality: 

oo 

dsP,7"$.(/) = ^A'")(r,t), (1.7) 

n = l 

where the monotone sequence of bounded positive operators A''"'(t, t) (see (1.4)) converges to if and only if the 
integral in the left-hand side is a densely defined operator. To make rigorous the above algebraic considerations of 
integrals, we must impose additional assumptions on domains and continuity. 

In the sequel we assume that the CP-map $*(•) is such that for each t £ R+ the map At{-) = A~^''^$t(-)A~^''^ 
is bounded and normal. In fact, the boundedness follows from the inequality A > $*(/). The Kraus theorem 
implies that any normal bounded CP-map At{X) on B{TC) {TC is a separable Hilbert space) can be represented as 
the sum A{X) = J2k -^k{t)XAk{t), J2k ■^kit)Akit) € B{H). This ensures a canonical representation of unbounded 
CP-map $4-) g: 

MX) = ^lii)X^k{t), ^k{t) = Akm\'\, 

k 

where A1.{t)Ak{t) £ B{7i). To study conditions sufficient for the minimal solution to be unital, we must extend 
the domain and continuity assumptions. We assume that for some N > 2 the operators A^^^^kit) are densely 
defined, 

Ojv C domAi'"'-l'fc(t), (1.8) 

and A^^'-^<i>k{s)Ws,ti' G L2°'=(R+,H) in variable s and norm-continuous in t. Thus the inequality (1.7) justifies the 
following assertion. 

Theorem 1.1. Under the domain and continuity assumptions, if the domain of the operator ds P^l'^9s{I) is dense 
in Ti, then the minimal solution of the Markov master equation (1.1) is unital. 

Since the sequence A'-"^(t, i) is positive and decreases monotonically, the sum 

oo 

C = ^ a„A'"'(r, t), a„ > 0, ^ a„ = oo 

n — 1 n 

converges to a densely defined operator only if A'"'(t, i) converges to 0. The series which correspond to this sum 
with On ~ can be represented as an integral of the minimal solutions of MMEs with the generators regularized 
as in 1^: 

Ct,x{B) = A$t(B) - GtX - XGt, A G (0, 1]. 
More precise, the series, representing the minimal solution of the equation 

^Pi,^'(B) = £,,.(P<^)(B)), Pi,^'(B)U=t = B 
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is the following; 



n = l •'^ ■'^ 



This identity and definition (1.5) imply 

nl ft 



f . . 1 

1^ d\ J dsPi,^)$4/)=X^iA(")(r,i). 



Therefore, the following assertion is true. 
Theorem 1.2. Assume that the domain and continuity conditions are fulfilled. If the operator 

"1 pt 



c = j^^ d\ j dsPiy^,{i) 

is densely defined in Tl, the minimal solution of the Markov master equation (1.1) is unital. 
Let us consider a priori bounds for the operator Ps,t^s{I)- 



(1.9) 



2 A priori bounds 

Assume that there exists a smooth family of positive self-adjoint operators At > $t(/) > / and N > 2 such that 
for any ip £T^n ^ domAt and 

$.,t(AO[V'] - 2Ke(GtV, AtV) - {^Atxjj) < ct||A,'/Vll, (2.1) 

(cf. Q and [^), where ct G I/i°'^(R+), Ct > 0. Such operator family is called a family of reference operators. We 
assume that the family of operators A~'^''^AtA~^/^, with the previously defined operator A, is densely defined on 
Ti and admits a continuation on the whole space TL which is uniformly bounded and strongly continuous. Let us 
prove that condition (2.1) ensures the a priori estimate 

^"/"(A^) < At e^- (2.2) 
This estimate can easily be proved by induction. Indeed, for all < r < t we have 

or OT ' 

> W:4^r{Ar)+CrArWr,t > M/,* ^ A, W^,, t = CrVr,tAr, 

Hence, by solving this terminal differential inequality, we obtain 

P^°^iAr) = V.,tA.<Atef^^'^\ 

Assume that P^"'(Ar) < AtcS* and let us prove this inequality for P^"+^'(Ar). 
From the recurrent definition of P^"^^-'(At) and assumption (2.1) we have 

Pi"t'^'\A.) = w:^tArWr,t + j'^ dsVs,t'i>sPi^"J{A^) 

< w:,tArWr,t + dse'- ''"'"'■Vs^t^siAs) 

< W:^tKWr.t + J' dseJ'' ""-"'-VsAcsA, + GIA, + A,G, + A,) 
= K.tA, + /* '^^ ^ (^■'^' "■''■K,tA,j = Ate^" "'^^ 

This estimate readily implies that the operator P™"(Ar) is densely defined, domA C domP™"(Ar), and hence 
the unital property holds. Therefore, the following assertion holds true. 

Theorem 2.1. Let the domain and continuity assumptions be fulfilled and there exist a reference family At- Then the 
minimal solution of Eq. (1.1) is unital. 
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Assume that for a formal generator £(•) with constant coefficients there exists some constant reference operator 
A, i.e. 

£(A)<cA, C{B) = ^{B)-G'B-BG, 

whore G — ^$(/) +iH, H — Hs + Hs.a., and let Hint ~ Hs.a. + A be a self-adjoiut operator. Then inequality (2.1) 
holds for Ct{B) = -GIB- BGt and the reference family At = Utku't, where Ut = e"^™* 

Gt = f + i{Hs - A)] Ut, MB) = u:^UtBU:)Ut. 



2 

int 



Let us discuss an opportunity to use some fixed reference operator A for problems with time-dependent coefiicients, 
which arise in the interaction representation. 

Let A > and Hs.a. be self-adjoint operators such that the sum Hint = Hs.a. -h A is positive and self-adjoint, and 
there exist fi, v > and e € (0, 1) such that 

-nHf^t < Hs.a < uHLt, < Hi„t < vk. (2.3) 

Note that for any positive self-adjoint operator X and e e (0, 1], we have Jf^ < 1 -\- X, since < 1 + A for any 
positive A, and hence 

X" = y^ \^Ex{d\)< J {l + X)Ex{d\), 
where Ex{dX) is the spectral family of the operator X. Then 

At = UtAU: = Ut{Hint - Hs.a.)U: 

= UtH-^/'Hs.aH-^/'u:^Hli 

< Hl^[l + I^H~^^"HLtHr^^'^Hl^ 

< Hint + H{Hrnt + /)<(!+ /i)j^A + //i < CoA. 

Thus under the above assumptions (2.3) there exists a constant co = p + (1 -h ii)v such that UtAUt < cqA. 
Assume that A is a reference operator for some formal generator £■{■). Then we have $(A) — G*A — AG < cA. 
Assume that 

HHs.a.) < ciA, ^{Hlt) < C2A, (2.4) 

and consider an estimate the action of the formal generator £(•) in the interaction representation o n the element 
A: 

£t{A) = u; (^^iUtAUt*) - G*A - AG^ Ut 

= u: (^^{Ut{Hint - Hs.a.)U:) - G*A - AG^ Ut 

= U: (^Hs.a. - UtHs.a.Ut*) + 4>(A) - G*A - Ag] Ut 



< U: \^<i{Hs.a. - UtHs.a.Ut*) + cAj Ut 

< u: (^^{Hs.a) + tiHHLt) + cA^ Ut 

< (ci + /^C2 + c)UtAUt < co(ci -|- /^C2) + c)A. 
Thus, under the above assumptions 

£t(A) < AA, A = Co (c-l-ci -I-MC2). 

By Theorem 2.1, this estimate implies that the formal generator Ct{X) — Ut jC{UtXUt)Ut is regular. On the 
other hand, the minimal quantum dynamical semigroup generated by jCt{-) is unitary equivalent to the minimal 
dynamical semigroup generated by _ 

£.{■) = £.{-)+i[Hint,-]. 

Hence £■{■) also generates a unital minimal dynamical semigroup, and it is regular too. Its coefiicients are $(•) 
(the same CP-map), and H = Hs — A + Hint = Hs + Hs.a.- Thus one can add a self-adjoint operator Hs.a. to any 
regular generator £.{■) which possesses a reference operator A if conditions (2.3)-(2.4) are fulfilled. In this case we 
call (A, Hint) a A-pair for the generator £{■)■ 

Theorem 2.2. Assume that the domain and continuity conditions are fulfilled. If for a formal generator £■{■) there exists 
a A-pair, the generator £(•) is regular. 
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3 Structure of generators of MME in quantum optics 

The typical formal generator £(■) of a Markov master equation in quantum optics (see (Schack & Brun '96), 
lH (Brun & Gisin '96), @ (Ariano & Sacchi '97), (ZoUer & Gardiner '97)) acts in B{H), H = (£2)®'^ ® C"; 
its Lindbladian form reads as follows: C{B) — $(-B) — CB — EG, where 



1 

$(B)=^$fc(B), G=-^I) + iH, 



2 

fc=i 

$fc(B) = Afea^Bafe or ^k{B) = ^jikUkBal, (3.1) 

As,, /Xfe > are positive operators in C*^, Ofc and a], are adjoint creation and annihilation operators acting on fc-th 
factor of the tensor product {t2)®'^ , i.e. ([ofc,a„] = 0, [a„,a^] = 5k,n), and 

H=Y. {^^ n ("D"^' («^)'"^' + (3-2) 
j ^ fe=i ^ 

is an operator in Ti represented by a symmetric polynomial of a finite degree in creation and annihilation operators 
with matrix coefficients hj G C*^ ® C*^ (see (Wiseman & Vaccaro '98), (Kist, Orszag, Brun & Davidovich 
'99). 

Note that the Hermitian structure (3.2) of the operator H does not imply its self-adjointness. For example, the 
Hamiltonian of the third order 

H = ((1 + i(« + - + (« - + i(« + a')'; 

is not a s.a. operator in I2, because it is unitarily equivalent to the symmetric operator H — i{{l+x^)dx+dx{l+x^)) 
in £2(9.), which has the nontrivial eigenvector 

Vl + a;^ 

such that Hip = — ii/;. Hence the symmetric operators H and H have the same nontrivial deficiency index, and 
TT = is a projector to the deficiency subspace X = Tid- 

Consider conditions on the projector tt which ensure the violation of the unital property for equations with constant 
operator coefficients. We recall that the condition necessary and sufficient for the minimal solution to be unital is 
the weak convergence to of the monotone sequence of bounded positive operators 

Q?(J)^0, Q,{X)^=^ f°° dte-''Vo,MX) (3.4) 
Jo 

(see Q). Hence the existence of a positive bounded operator X, \\X\\ < 1, such that Q^iX) > X for some e > is 
sufficient for the violation of the unital property, since the sequence 

Q"il) > Qe{X) >X>Q (3.5) 

clearly does not converge to 0. 

Theorem 3.1. // the coefficients of a formal generator C{-) satisfy domain and continuity assumptions and there exist 
a positive bounded operator X, \\X\\ < 1, and e > such that 

£.4X)[i>] > eX4i>] VV- G domG^ = Oiv, (3.6) 

then the corresponding minimal quantum dynamical semigroup does not preserve the unit operator. 

Proof. Let us derive the inequality Qe{X) > X from (3.6). Inequality (3.6) implies that ${X), > {eX + G* X + 
XG)t on T>N- Then for ij) G T>n, we have ipt = Wt?/' G Ojv and 

e-"*-l>.(X)[Vt] > e-"*(e(Vt,X^/>0 + (GV't,Xi/>t) + {XiPt,GiPt)) 

^-i{e~-\\XM?). 



Therefore, 

' e-^"$(X), [V-r] dr > ii:, XiP) - e-'\iPuXiPt). 



f 

Jo 
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The limit as t ^ oo yields an inequality for X: {tl),Qe{X)tp) > {tp,Xtp) for any ^ £ Dn by definition of the map 
Qe(-). Since Dn is dense in Ti. and the map (5e(-) is bounded, this inequality is equivalent to 

Q.{X) > X, (3.7) 

which contradicts the necessary unitality condition. □ 

A natural candidate to be used as X in inequality (3.6) is the projector to the deficiency subspace of the operator 
H, if such a subspace exists. 

Theorem 3.2. jl^ (Chebotarev & Shustikov '00) 

Let H be a densely defined symmetric operator. Assume that it has a nontrivial deficiency subspace TCd = {i^ ■ H*i' = 
—i^} and TVd is the projection onto T~Ld- If moreover, there exists e > such that 

<E>(7rd).[V'] -5Re($(/)iA,7rd7/)) > -(2 - e)||7rdV||' (3.8) 
for all ip G I?jv, then inequality (3.6) holds and the necessary unitality condition (3.4) is violated. 
Proof. Let us prove that (3.8) implies (3.6) in the following form: 

.C(7rd),M > elkdVIl' VV'G^iv. 
Indeed, since H*TVd = —ii^d and tt^ = ttJ, it follows from (3.8) that 

C{7Td)M = $(7rd),[V'] -5Re mi)i,,7Vdi^)+i{{H'4^,7Vdip) - {-Kd^^H^P)) 

= HTTd)*\4}] - (<E>(/)V',7rdV) + 2||7rdV||' > elkd^f (3.9) 

and inequality (3.6) is true. □ 

Note that for any finite polynomial H — H2 + Ha. a. in creation and annihilation operators a], and Ok, there exists 
a diagonal operator 

A = CA('l + ^(aIafc)'"'=Y CA>0 (3.10) 

such that H2 and Hs.a are relatively bounded by A with the relative upper bound 0(c^^). One can use A as the 
reference operator. In any case we assume that $(/) > / is a s.a. operator and 

adomA C domA^^^, dom A C dom A^''^. 

The last two assumptions readily hold if mk > M — 2 in (3.10). 

Theorem 3.3. // the Hamiltonian H can be represented as H = H2 + Hs.a., where Hs.a. is a self-adjoint polynomial of 
a finite order M in creation and annihilation operators and H2 = H2{a\a) is a polynomial of the second order, then 
there exist ca > and {mk} > M such that (A, Hint = Hs.a. + A) is a A-pair for the generator (3.1)-(3.2). 

Proof. We recall that any finite polynomial in creation and annihilation operators of order AI can be dominated by 
the diagonal operator (3.10) of higher order, provided the constant cd is sufficiently large and TV = min{mfe} > M. 
Hence for sufficiently large c_d, and A*' by the classical perturbation theory |l^]. Hint ~ Hs.a. + A is a positive s.a. 
operator such that dom Hint ~ domA, and Go ~ i{H2 — A) + $(/)/2 is an accreative operator, domGo = domA. 
Since <!>(•) and Hs.a. are operators of a finite (second) order, the property (2.3) and (2.4) of A-pair can readily be 
fulfilled by choosing M, N and ca sufficiently large. 

The commutator of a polynomial of the second order in creation and annihilation operators with arbitrary poly- 
nomial of order A4 < 00 has the order M or less. Hence, the commutator 

i[Ho,A]^i[H2-A,A]^i[H2,A] 

is an operator of the same order as A, and hence there exists a constant c £ R such that i[Ho,A] < cA. 

A simple algebra shows that for CP-map (3.1), the operator $(A) — (A<1?(/) -I- <I>(/)A)/2 is also a polynomial of the 

same order as A. More precise, the following two estimates hold: 



{alYAai ^ i (^{alYaiA + A(4)'aL j < 0, 



for any I > 0, and on the other hand there exists c > such that 

OkAal - {okalA + Aokal) /2 < cA, 

for the operator A (3.10). Therefore, there exists c > such that £o(A) < cA on domA, and hence A is a reference 
operator for 

$fc,i(B) = Xk.iialYBai, $^(5) = XkakBal 

This proves the theorem. □ 
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Examples 

In this section we consider some classes of Hamiltonians and completely positive maps for which our Theorem 2.2 
is applicable. 

1. Let A be a complex number and m, n > 0. Set 

f/ = A(a|)'"a^ + Aar(4)", A G C. (4.1) 

Let us prove that all Hamiltonians of such form are essentially self-adjoint in Ti.2 = ?2 <8) ^2. It suffices to prove that 
there does not exist a vector 

iP = {V'fc.j, k,j > 1, ^ \ipk,jf = \ W\n2} 

k.j 

such that H'tp = ±ii(j We set t/j^j = if min{fc, j} < 0. 
Let us rewrite these equations for components iptj as follows: 

±ii,k,j^XA^fi,k-,r.,j+n+XBiyi,k + n^,,-n, A^f , B^f > 0. (4.2) 

We set tpkj = if min{fc, j} < and skip exact expressions for the functions A^'J^ and B^'" because they are 
irrelevant for the proof. The important property of this system is that it splits into a set of independent Enite 
subsystems of linear algebraic equations with respect to values of one of the components of the set 

Xk = {xj = i>k-jm,i+jn, k-jm> 1, jn > 1, j = 0, 1, . . . , [k/m] - 1}, 

where xj = for all j < 0. For each k, m, n fixed, the system of linear algebraic equations corresponding to (4.2) 
has the three-diagonal form 

±ia;j = \AjXj+i + XBjXj-i 
with some positive Aj and Bj. But it is a well-known fact (see jl^]) that 

/ ±il XAi ... \ 



det 



A*Bi ±il XA2 



(4.3) 



\ ... X'Bn-i ±il I 

where the entries of the matrix are ik x fc)-blocks, A and A* are Hermitian adjoint (k x A:)-matrices, and I is the 
unit matrix in C*^. 

By the Gershgorin theorem Hamiltonians (4.1) are relatively bounded by the diagonal matrix A (2.3) of order 
M > m -I- n, and the relative upper bound decreases as cu — ^ 00. Hence all formal generators with the completely 
positive parts (3.1) and Hamiltonian part (4.1) are regular. 

2. The same assertion is true for Hamiltonians from C{lf^ ^ C'^') of the following form: H — Hint + ^^o, 

H,r.t = X{a\r'aJl' ■ ■ ■ a™"(ajv)"" + X'a^\a\r' ■ ■ ■ (a]v)'""Q,;" , (4.4) 

where A and A* are Hermitian adjoint (M x Af)-matrices, ^ nrik + Uk — K, and Ho is any symmetric operator 
dominated by A and such that 

3ceR: i[Ho,A]<cA. 

The proof of self-adjointness of Hint is based on a similar factorization of the set of block-matrices {i/'fei,...fcjv} ^ 
if^ ® C*^ and on the reduction of the homogeneous system of linear algebraic equations to the set of finite- 
dimensional linear equations with nondegenerate three-diagonal (Af x M)-block matrix (4.3). 
As in the previous case, the interaction representation is generated by the self-adjoint operator Hint dominated 
by the diagonal operator D for M > X]('^fc + ^k), and the Hamiltonian Hq of £o{-), because it satisfies the 
conservativity and compatibility conditions. For generators (1.1), any symmetric operator on Ho of the second 
order in creation and annihilation satisfies the above assumptions. 

3. Consider the physical example (Schack, Brun & Pecival '96) of a formal generator £{■) in B{l2 (8) Z2 ® C^) 
with CP-part (1.1) and the Hamiltonian 

H — Ei{al — ai) + ^i{a\^a2 — aial) + ua+a- + r]i{a2(J+ — a\(T-), (4-5) 

where E is the strength of an external pump field, x is the strength of the interaction, uj is the detuning between 
the frequency of the field mode 02 and the spin transition frequency, and 77 is the strength of the coupling of the 
spin to the field mode 02. The completely positive part of the generator reads as follows 

$(B) = 2'yia\Bai + 2'y2alBa2 + 2na+Ba.. (4.6) 
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It describes the dissipation of the field modes and the spin with coeflicients 71, 72, and k, respectively; cr± are two 
by two matrices. 

The Hamiltonian (4.5) can be readily represented in the form (4.3) with k = 2, K — 3, X — I, Hint ~ ^i{cil^ci2 — 
fliOj) and Ho — H ~ Hint- The completely positive part has the form (3.1). Hence the formal generator (4.5)-(4.6) 
is regular. 

4. The kinetic stage of the evolution of a quantum system interacting with environment is described in (Kilin 
& Schreiber '97) by the following Markov master equation: 

^ = -iu \H{a\a), 

+ T2{n2 + 1) { [aV, (at)^] + [{a^f, aa^] } 
+ V2n2[[{a^fa,^'\ + [a^a(at)2]}, 

where s.a. operator H — H{a\ a) is a finite symmetric polynomial in and a of order no greater 4, r2 — TTK^g2 
is the decay rate of the vibrational amplitude. Here, the number of quanta in the bath mode n2 = n{2ui), the 
coupling function K = K{2uj), and the density of bath states 32 = g{2uj) axe evaluated at the double frequency of 
the selected oscillator. The corresponding dual CP-map <3>(-) acts as follows 

t^2 EJ„2 , ^ 2 D/-„t\2 



$(B) = 2r2 (^(n2 + 2){a^yBa + n2a B{a 

This case is rather simple: A — c(at)^a^ and Hint = H + A, where c is sufficiently large: A > "!'(/), and \ \Ah\\ > 
2\\Hh\\, so that Hint is a s.a. operator, provided H is s.a. operator. Hence the generator of the above master 
equation is regular. 

5. The previous example can be generalized as follows. Set 

$+{B) = {a^rBa^, $" (B) = a™B(at)™, A„ = ((a+a)" + /)A, 
A > 0. Then we set ^ 

An^N = A(iV" + 1)^^, <^iiI)^N = ( ^^ir^- ) ^ 



V iV! 

for TV-particle component of i/^jv of the vector tp — {i/jq, V'l) . . . } G '2. Hence there exists n > m and A — A(m, n) > 
such that An > $m(-^)- 
Similarly we obtain 



<1>±(A„) - (A„<1>±(/) + $±(/)A„)/2j^^ = 

= A( — 1 ((iV=Fm) -N )i:N. 

Therefore, for any formal generator with completely positive part 

$(B) = ^ [B) + c^$-, (B)} 

k 

with positive matrix coefficients € C*^ ® C*^, the third A-pair assumption is fulfilled if the balance condition is 
true: 

-pE{4(^^[(i-Wiv)"-i] + 

+'^-.^-^^[i^+n./Nr^i]}<ci. 

Then the regularity conditions of Theorem 2.2 are fulfilled for A sufficiently large if H = H2 + Hs.a., where H2{a, a^) 
is any Hermitian quadratic polynomial, and Hs.a.i^j'^'') is any self-adjoint Hamiltonian of order less or equal 2n. 
The balance condition is readily fulfilled if 

S = ^{cjmfc - c^nfe} > 
fe 

is a strictly positive operator in C*^ (g) C^'^. 
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The generators of MME with completely positive component of the fourth order was used in (Schneider & 
Milburn '97): 

-H = ^2®C^ 'i>{B) = {a^a+ + {a''fa^)B{aa^ + {a''fa+), (4.7) 

where 



Hence, for diagonal operators B £ B[TL) 





^ 1 \ 




' 






, j 






I) 



B=[ I' ), B,,,eB{h 



we have (t± = and 



a^Bi{a'^f 
{a'^fB2a 
Therefore, the generator C[-) has the component 

£ii(X) =a'X(at)2 - {a^ {a'^f X + X{a)f) /2 + t[H-,^, X] 

which is unregular for any first order operator i/n. 

6. The paper (Lanz, Melsheimer & Vaccini '97) presents examples of formal generators inli. = h^h with the 
CP-component $(-B) =0103-8 02 a|. Let us prove that for the generators with coefficients 

^(B) =af(at)"Baf (at)^, L,M>1, H = H2 + Hs.a.{a\a) (4.8) 

there exists a A-pair. Consider the generator 

Co[B) = ai [a!2) Baj {a[) -Boai{a[) (a^) . 

Straightforward computation proves that $(/) is not a reference operator for -Co(-)- 

Lemma 4.1. For any N > 0, there exists a polynomial 

N 

A^ = AWj + 5]Ai^'(al)^-Sf-'=(at)%\ A(^\ a'^' > 1 (4.9) 

fc=0 

and a real constant cn such that /lo(Ajv) < cjvAjv. 
Proof. Note that 

£o{(a[) 01(02) a2)=mL{a[) (aJ,) 

— nM{a[) a-i (ttj) 02 -|-l.o.t., 

where the lower order terms (l.o.t.) can be dominated by the main terms and A'^'/ for all A*-^-* sufficiently large. 
Hence, 

JV-l 

C,{Kn) = [L{N - fc)A(^' - M{k + 

X (a|) (aj) + l.o.t. 

Therefore, all main terms have negative coefficients if 

,(N) ^ ,(jv) [N -k)L 
{k + 1)M 

For the fixed A*', the lower order terms can be dominated by the main terms plus A''^'/ for all A'-^' sufficiently 
large, that is 

3c>0: Z:o(Ajv) < cA'^'/ < cAjv 

if (4.10) holds. Since 

A \(N)r I UN N , \ (JV) / UN N 

where the coefficients A'^', Aq^\ X\^^ can be chosen greater than any constant c > 0. In particular, they can 
be chosen such that the diagonal s.a. operator Ajv dominates with arbitrary small upper relative bound a given 
polynomial Hs.a.- Note that for any quadratic operator H2 = H2{a\ a), Co{H2) is a symmetric polynomial of order 
2(1/ + M) in creation and annihilation operators. Hence it can be dominated by Ajv for any TV > 2{L + M). In this 
case, {Ajv, Hint = Ajv + Hs.a.} is a A-pair for the formal generator £(■) with coefficients (4.8). This proves that 
£(•) is regular. □ 



^iz>>^r'ir-T^- (4-10) 
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5 Discussion 

By using the concept of A-pair, we have analyzed the regularity property for a wide class of generators of MME in 
quantum optics which are available for authors. We proved that the generators of the form 

C{B) = - ($(J)B + B$(/))/2 + i[H2 + Hs.a.,B] 

are regular for CP-maps (3.1), (4.6), and (4.8) if Hs.a. is a self-adjoint polynomial of a finite order in creation and 
annihilation operators, and H2 is a symmetric operator of the second order. To conclude the paper, we recall the 
most important open problems. 

Generator of MME can be irregular if $(■) is as in (4.7). From mathematical viewpoint, to select a unique solution, 
one must introduce a kind of boundary condition as was done in |2^], where all unital extensions of the minimal 
quantum dynamical semigroup are described in terms of extension of its resolvent. In analogous classical cases, 
the boundary conditions for stochastic processes follow from Dynkin's formula for infinitesimal operator 

of the Markov semigroup. The physical sense of boundary conditions for quantum systems should be related to 
conservation laws, but physical examples of MME with boundary conditions still are not known. 
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